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ABSTRACT 


This paper investigates the feasibility of an air-to- 
surface missile delivery system employing impact prediction 
to assist a human operator in controlling the system. The 
impact prediction information was used to drive a display 
which the operator used in controlling the flight of the 
weapon. It was desired to find whether impact prediction 


proved useful in weapons control. 
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I. INTRODUCTION AND STATEMENT OF THE PROBLEM 


A. TYPICAL ASM WEAPONS SYSTEMS 

Most older air-to-surface weapons share the characteris- 
tic that once the weapon is launched it is under the influence 
of only the ®environment, e€.g., gravity, wind drive, saraemere., 
Examples in this category of weapons are bullets, ageneral- 
purpose bombs, and rockets. This makes it the responsibility 
of the weapons officer to insure that the fire-control solu- 
tion is accurate enough for the weapon to impact at the tar- 
get within the kill radius of the warhead. In the normal 
tactical situation this is a difficult task to perform prop- 
erly. As a result some of the newer weapons systems have 
been designed with a capability for guidance of the weapon 
to the target. 

These newer systems require that to provide guidance for 
the weapon, the weapons control officer must maintain visual 
contact with the weapon during its flight. Again, this solu- 
tion to the problem suffers from tactical restrictions. In 
a combat environment it is very difficult to maintain a flight 
path which allows the weapons control officer to maintain 
visual contact with the weapon during its flight. 

In order to overcome the shortcomings of these systems, 
even newer systems such as the Phoenix employ very sophisti- 
cated missiles which allow the launching aircraft to stand-off 
from the target and still deliver a weapon so as to destroy 


Or seriously damage the target. One aspect of these systems 


which is objectionable is the high cost of the missiles which 
are employed. As a result this study was made of a system 
which would employ a stand-off delivery capability while em- 
ploying a relatively unsophisticated, and therefore inexpen- 


sive missile to deliver the warhead. 


B. ANALYSIS OF STAND-OFF WEAPONS DELIVERY 

Stand-off weapons delivery is best defined as the ability 
to remain outside the destructive range of the enemy's weapons 
and deliver ordnance on the enemy installation. The advan- 
tages of remaining outside the destructive envelope of the 
enemy'S weapons are obvious. The disadvantages of standing 
off from the target are numerous. The greater the stand-off 
distance, the more accurate the guidance and system signals 
must be as well as the knowledge of the target location. 

This is necessary to provide a reasonable probability of tar- 
get kill. 

In its most basic form the air-to-surface weapons delivery 
problem consists of two parts: the detection and location of 
the target by some means, and delivery of a weapon to the 
target. 

The solution which is proposed in this paper is based on 
three assumptions about what the author feels an air-to- 
Surface weapons system should incorporate. These three assump- 
tions are 1) that the target should be located by some method 
which does not endanger the launching aircraft; 2) that the 
primary purpose of air-to-surface missiles is to deliver ord- 


nance on the target, not to carry a large amount of expensive 


electronic equipment to be blown to pieces along with the 
target; and 3) that to improve the overall accuracy of the 
system, the missile will be under some form of command guid- 
ance during its flight to the target. 

The solution which is proposed in this paper is one of 
the "shoot and loiter" school of weapons control. Those of 
this line of thought contend that it is more economical to 
shoot a relatively unintelligent missile at the target and 
maintain the sophisticated system parts in the launching air- 
craft. This requires that the launching aircraft remain in 
the target area until the weapon has hit the target which 
increases the probability of the launching aircraft being 
destroyed. 

The opposite point of view is that the aircraft should 
launch the missile and get out of the target area. This 
approach necessitates the use of a very sophisticated missile 
but reduces the probability of the launching vehicle being 
eesStroyed. 

The trade-off to be considered between the two approaches 
is to obtain the proper balance of economy versus probability 


of survival of the launching aircraft. 


Ii, THE, PROBLEM: ANALYSTS AND (SOLUMION 


In) this chapter the two divisions of the problem willie 
discussed in detail. The first is the problem of target loca- 
tion. The second, delivery and guidance of the weapon, will 
be subdivided into two sections. The first section will be 
concerned with the initial fire-control solution for the mis- 
Sile flight path and the second with the guidance of the 


missile to the target. 


A. LOCATING THE TARGET 

The task of locating the target will be limited to the 
case of a radiating target such as a radar. This limitation 
was made in order to propose a solution in which the launch- 
ing aircraft both locates and destroys the target. The 
assumption of a radiating target allows the use of a direction- 
finding method to locate the target. 

Several methods of direction finding are available. These 
include triangulation of radiation given off by the target 
using two or more receivers at known locations, the use of 
filtered sequential measurements of bearing and depression 
angles from a single aircraft at known. positions, or even. tie 
Simple method of having two ground observers locate the tar- 
get. For this solution the second method was used. 

Poot ome oOmeury gon langen loOcacon 

In Fig. 1 the basic principle of location by triangu- 


lation is illustrated. The points 1 and 2 could be either 
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Figure 2. Definitions of Azimuth and Depression Angles 
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por cdis ClNnCEmaccemvers lis teningirona signal at the same time, 
or they could be the two locations of a single aircraft at 
different times. 

Figure 2 gives a definition of the angles involved in 
all the target-location schemes. The azimuth angle is that 
angle measured from a true north reference positive clockwise 
to the line-of-sight (LOS) between the aircraft and the target. 
The depression angle is the angle measured from a reference 
plane normal to a line connecting the aircraft and the center 
of the earth. An angle which rotates below the normal plane 
is considered negative and one which rotates above the plane 
is considered positive. Appendix E contains a discussion of 
the reference systems used in the entire problem. 

In the following two sections two different methods 
of location of the target by use of angular information will 
be presented. 

Zo caeron Using POsttiOnalgand Angulargintoemation 

The first of the two methods for target location will 
be explored here for the two-dimensional case to show the 
principle involved. An analysis of the three-dimensional 
case is contained in Appendix A. 

Figure 3 shows the geometry of the two-position method 
of target location. The aircraft, at two different times and 
at different positions, measures the azimuth angle from the 
aircraft to the target. At the same time as the angle is 


measured, the position of the aircraft 1S recorded. Assume 
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that the target coordinates are xy t 


are listed below. 


Time Aircraft Location Azimuth Angle 
1 Xy Yy Oy 
Z Xo Y5 0. 


From Figure 3 it can be seen that 





eK 

Tan O4 = roy (TI-1) 
eel 
ee 

Tan 9,5 = <3 (1 Eee 
tee 


Writing the two equations as simultaneous equations 
and solving them: 


Tan 90 


ce 1 it 
(II-3) 
xy = x. = (Y,-Y.) Tan oe 
x, - vy Tan oF = xy = vy Tan 04 
(II-4) 
x ~ Vy Tan 04 = x. = Y~ Tan 0. 
Letting A = Tan Oo, B = Tan Oo, C = xy _ Yy Tan O17 


and D = X, - Y, Tan 6, and writing the equations in matrix 
form 
iee= X, E 
= (II-5) 


Using normal matrix algebra technigues, the solution 
can be shown to be 


(II-6) 





e A-B : (II-7) 


Expanding the results into their fue romieby seoclacing 
the A, B, C, and D with their respective values gives 


Tang -Y.Tané 


1% 5 >)-Tané 


> (X,-¥,Tan6,) 


ee ee eee (IIT-8) 
t Tan, Tang. 
and 
~ (X,-Y,Tan68,)-(X,-Y,Tané, ) aa 
i Tan9,-Tan6., ° 


This technique of locating the target can be quite 
effective and accurate if the two positions used in deriving 
the solution are far enough apart and if the bearing angle 
measurements are respectably accurate. Ideally the two bear- 
ing angles should be separated by ninety degrees. This sep- 
aration will reduce to a minimum the effects of any noise in 
the measurement of the angles. If the target radiates for 
only short periods of time and the ambient conditions are 
such that the measurements are very noisy, the derived target 
location could be greatly in error. For this reason an al- 
ternate approach to the problem of locating the target was 
taken. The alternate method uses only angular information 
and a hybrid filtering technique. 

3. Location Using Angular Information and Hybrid Filtering 

The target location method discussed in this section 
is the result of research done by Demetry [1]. The method em- 
ploys an unusual hybrid filtering technique in which a Kalman 
filter is considered linear for calculation of a gain sched- 
ule but in which prediction is done nonlinearly. The hybrid 


filter will be discussed here briefly. 


tS 


It is desired to locate a radiating target by a method 
in which only a bearing angle to the target is measurable. 
The equations below give the relationships used in the formu- 
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el 
a= Rk (II-11) 
trek 
is the estimated target location at time k and 
x(t) 
Bey) 
= jy(t) gl igd il) 
aise 


e. 
is the aircraft vector. 
The aircraft position and velocity components are 
considered to be measurable without noise for this evalua- 
tion as they are throughout the remainder of the problem. 
The bearing angle dynamics are described by the 


differential equations 


Oat) 
w(t) = [6(t)| = £(a,&(t) ,t) (II-13) 


6 (t) 
PieGoe i eicoeadsiOn near vector function. The angle measurememe: 
at any discrete time k is 


ar 


HY, + v,, k = = (aed) 


ipemen JH =i) “O50 ) and Vv, is the error in the bearing angle. 
The Vy 1s random in nature with an assumed mean value of zero. 


The gain schedule was calculated by assuming linear 


dynamics and constant acceleration for the process and using 
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the recurSive Kalman gain equations [2]. However, the predic- 
tion of angular information was done by using the actual angu- 
lar dynamics in a discretized form rather than employing the 
linear ® matrix which was used in computing the filter gains. 
These equations, by using true system dynamics to make approx- 
imate predictions, reduce the location error since, in effect, 
they consider the system linear. In this system this is not 

a bad assumption since the sample period is relatively short 
compared to the rates of change of the system variables. 
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This is a brief discussion of this particular method 


of location of a target. For the reader interested in pursu- 


ing this further, Ref. 1 will be of interest. 


Peli PtRE-CONTROL SOLUELON 

Before it is possible to launch a missile at a target it 
is necessary to develop a solution for a flight path which 
will take the missile from the launch point to the target. 
The accuracy with which the solution need be computed is en- 
tirely a function of the amount of guidance which will be 


available to influence the missile during its flight. Also, 


te 


before an initial solution can be computed, it 1s necessary 
to know what approximate flight path the missile will follow 
from the launch point to the target. 

In the solution which is being investigated here, the 
proposed flight path for the missile is a nominal parabolic 
path. The path the missile actually travels will not be para- 
bolic due to air drag, control effects, and other outside 
influences. 

There are two main reasons for selecting a parabola as a 
nominal flight path. The first concerns a relationship be- 
tween the height of a parabola and distance to the target 
from a point under the peak of the parabola. This ratio is 
referred to as the A/R ratio and its importance is discussed 
in Appendix B. The second consideration is that the missile 
will be under the influence of gravity in a very helpful man- 
ner. The gravitational force will help to maintain the veloc- 
ity of the missile and hence, the controllability of the 
missile, after the boosted phase of the missile's flight. 

ie CCOMC EL VAC meno gh) ne Controlsco lution 

The general geometry of the fire control problem is 
illustrated in Fig. 4. In the solution the aircraft, point 
A/C on the figure, is considered to be at the X origin of the 
figure and the target, point T in the figure, is considered to 
lie on the X-axis. The parabola will be defined such that it 
contains the two points of the launching aircraft (A/C) and 
the target (T). Since a parabola is defined by a second-order 
polynomial, e.g., of the form y = ax? - Dx + cy come “thing 


point on the parabola must be defined. 
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Figure 5. The Nominal Parabolic Flight Path with 
Distances Defined 


V9 


In Fig. 5, the parabola of Fig. 4 is redrawn, this 
time witn three points and four distances defined. For the 
purpose of defining the nominal parabola, the launching air- 
craft is placed on the Z axis so that range a zero at the 
Origin. The altitude of the aircraft is labeled ZA. A sec- 
ond point is the location of the target. The target is at 
zero elevation and a distance D from the origin. The third 
boume, wach 1S necessary for the parabola to ber detined 7 yr. 
the peak of the parabola. The altitude of this point is A 
and it is D - R from the origin. Using these three points 
it is possible to solve for the coefficients of the second- 
order polynomial defining a parabola. 

Writing the equations of the parabola at the three 


points gives 


a(0)* + b(0) +c = ZA (II-18) 
ip Re te (D Reick (II-19) 
a(D)“ + b(D) + «¢ #0 (II-20) 
R= kD : (ILI-21) 


Equation II-21 is used to define the relationship 
between the range from the origin to the target and the 
range from under the peak of the parabola to the target. 

2. Method of Obtaininggthe Launeheingle 

In a system in which the launch point and target point 
are coplanar in the earth plane, a ballistic mass will travel 
furthest when the velocity vector is initially directed 45 
degrees above the horizontal. This.is true as long as the 


mass has no external forces acting on it other than gravity. 
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Two factors involved in this investigation make the 
choice of a 45-degree launch angle unacceptable. The first 
is that the launching point and target are not earth coplanar 
and, secondly, that the mass is subject to thrust from a 
rocket motor. Thus, the problem of computing the launch angle 
becomes one of maximizing the range as a function of the dis- 
tance between the target and launching point in the Z plane. 

The results of this particular maximization were 
interesting. For aircraft locations which were uv to 5000 
meters noncoplanar with the target, the launch angle for max- 
imum distance was found to vary only one degree, between 47 
and 48 degrees. (The launch angle as used here is defined as 
the angle between the velocity vector and the horizontal at 
the launch point. This is not the angle of the velocity vec- 
tor at the time of burnout which is the angle which deter- 
mines the maximum distance of flight.) The launch angle was 
set at a constant 47.5 degrees as a result. This angle should 
give the maximum range regardless of range to the target. 
Starting with the launch angle that gives maximum range should 
prove entirely acceptable since one of the assumptions about 
the system was that the missile would be under some type of 


guidance during its flight 


C. GUIDANCE BY THE HUMAN OPERATOR 

Tnroughout this paper mention has been made of the fact 
that the missile will be guided by the human operator through- 
out its flight. In actuality the missile was not released to 


the operator's control until after it had finished the boosted 


Zul 


pitese Gf ats flight. During this time gehe system turned the 
missile up to the launch angle and yawed it to the prover 
flight neading. Once the missile started its free flight, 

a display of the predicted impact point saeies missile rela- 
tive to the target was generated for the missile control 
officer (MCO). At the same time, the control surfaces of the 
missile were set to a neutral position and control of the 
missile was given to the MCO who then guided the missile the 
remainder of the distance to the target. 

In order for the predicted missile impact point to be 
displayed, it was first necessary to track the missile and 
based on the tracking information, predict its impact point. 
The prediction information was then presented to the MCO in 
such a manner that he could make decisions as to the guidance 
control to be given to the missile. 

Peel roc agmesesnissile in Firone 

In this problem the missile was tracked in flight by 
a radar which was capable of providing four items of informa- 
tion. These were the range to the missile, the range rate of 
change to the missile, the azimuth angle to the missile, and 
tne depression angle to the missile. The reference system by 
whicn the angles were measured is discussed in Appendix E. 
The range was measured in the positive sense while the range 
rate was considered positive for a missile whose range was 
increasing and negative for a missile which was closing the 


iPenChindguamnoraft., 


Ze 


Since the radar was operating in a natural environ- 
ment, the information provided by the radar was not perfect. 
This noise would introduce errors in the estimation of the 
location of the missile during its flight, so it was desirable 
to use some technique to reduce the effect of the noise on 
the measurements. The technique used was Kalman filtering. 
Before deriving the necessary Kalman filter for this problem, 
the desired missile information vector should be investigated. 

Figure 6 presents the geometry for developing the 
missile information vector from the spherical vector of the 
radar. The desired missile information vector contains s1x 
components: the XYZ cartesian location of the missile rela- 
tive to the aircraft and the XYZ cartesian velocities, again 
relative to the aircraft. The acceptability of providing 
only a six-element state vector for the missile is dependent 
upon the sample rate at which the state vector is developed. 
Even though the missile is being acted upon by external forces, 
it is possible to consider the missile to be a piece-wise 
constant-velocity missile in either spherical or cartesian 
coordinates if the interval between samples is sufficiently 
small. 

From Fig. 6 the equations for the missile location 
relative to the aircraft can be developed. Then by the use 
of the time derivatives of the equations of location, the 
equations of the cartesian velocities can be found. First, 


the equations of location. 
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Zo-= Resim 6 (Liao) 


m 
x = R' sin 8 = R cos $6 sin 6 (ITI-23) 
= R' cos 8 = R cos 6 cos 86 (ITI-24) 


For the case shown in Fig. 6, the angle ¢ is negative, 
giving the proper negative value for Za Likewise, the signs 
of the remaining two values will be correct depending on the 
quadrant of the angle 0. 

To form the equations of the velocities in the carte- 
Sian frame, it is necessary only to take the time derivatives 


of equations I1I-22, I1I-23, and II-24. These are given below. 


Ze = R cos 6 r - R Sin @ (II-25) 

‘7 = R cos ¢ cos @ 6 - R sin ¢? sin @ ; + R cos ¢ sin 86 
(Ti=26y 

v = -R cos ? sin @ Q - R sin $6 cos @ 4 + R cos ¢ cos @ 
(II-27) 


From equations I-22) cheuga (il —2 7 ag tS eo meee 
necessary spherical values of missile movement to develop 
the cartesian state vector can be collected. These equations 
show that it 1s necessary to have a vector of spherical values 
of range, azimuth, depression, and the associated rates of 
change of the three quantities. From this requirement a Kal- 
man filter can be developed which will not only reduce the 
effect of the measurement noise on the quantities which are 
Peeing Measured;ebut also provide estimatesmet Ge eaneeeea. 
rates which are not measured. 

The Kalman filter for filtering and estimating the 


radar vector has the general Kalman predictor-corrector form 


2 


of equation II-28. The G, 1s computed from the recursive 


equations of I1I-29, II-30, and II-3l. More complete details 


@Qeeciie Cerivation of the filter and Ene vawewes used an the 


gain schedule computation are given in Appendix C. 
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Kel = Sefer + Se (BR ed) i 
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Baie = (I-GyH) Py yy (II-30) 
eae = OP. 1? + Q (II-31) 
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In predicting the missile impact point it is desirable 
to have an algorithm which is fairly accurate. This is the 
case Since it is the point at which the missile is predicted 
to strike that will be used by the MCO to develop his guidance 
commands. A requiewiGemniememsce si Gaeemeacy over a ballistic 
model can be obtained by introducing the effects of drag into 
the equations of motion for the missile. For the prediction 
algorithm developed here, the drag was considered to be pro- 
portional to the velocity of the missile. 

In this solution the drag coefficient of the missile 
is taken to be constant from the time of the prediction until 
the missile strikes the target. This, of course, is an invalid 
assumption but one which is necessary in order to develop a 
closed-form solution for the prediction equations. The error 
which is introduced by this assumption will have the effect 
of causing the MCO to hold the missile aloft longer than would 
normally be necessary. This action by the MCO will tend to 
increase the A/R ratio, thus making the missile impact on the 


26 


target more vertically and hence less subject to error from 
target altitude location error. (See Appendix B). 

a. Equations of Motion for Impact Prediction 

The equations of motion for impact prediction were 

written for each of the three component directions of the 
cartesian reference system. The equations were based on an 
assumed point-mass model of the missile subjected only to 
drag and gravity forces. The X and Y equations are similar 
to each other, varying Only in the drag coefficient. The Z 
equation takes into account the effect of gravity. The equa- 


muons Of mMmotLon are 


mk = -K,X CTresoy 
mY = -K5Y (Tae 3) 
mZ = -K 2 - mG (II-34) 
where 
K, = o|X|A, (II-35) 
K, = o|Y [AL (II-36) 
K, = o[ Za, (II-37) 
and m = missile mass, G = gravity, and A. = reference area. 


The method for predicting the impaece voine OF seme 
missile is to first compute the drag coefficients (Ky, Koy 
and K 2) which are then considered constant until impact. The 
Z motion equation is then solved for the time when the missile 
Z% location is equal to the target Z. Pirst; @teis neeeseae, 


to solve equation II-34 for the missile Z as a function of 


abies 
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b. Impact Prediction Equations 
The result of solving for the Z motion as a func- 


t10On OfStimewes 





ACE) =. (Beaver? = peel (Ilse 
where 
i 2 
Sa Zao)) m G 
3 K 
B= MG (Trae 
K 
3 
Ci — Z_ (0) (II-41]) 
icy | 
m 


This mixed equation is not directly solvable for 


tiawevernuie Of Lime Such that Zi 1s equal to Z The Newton- 


be 
Rapnson iterative method is employed to find the value of 
time which will satisfy the equation for a specific set of 
conditions (a particular Z 69) and z_(0)). The value of time 
thus found is substituted into the X and Y solutions to find 
the predicted impact point of the missile. The X and Y 


predicted positions at impact are given by equations II-43 


and II-44. 
mX_ (0) -K,t, /m 
i = XO) + K, (l-e ) (II-43) 
can (0) eh t./m 
Y = Y_(0) + ae (l-e ena ) (TI-44) 
p m K. 


where t. is the value of impact time from the solution of 


equceion Iil-33., 
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In order for the MCO to make intelligent decisions 


regarding the application of control to the missile during 


flight, there had to be some method for displaying the pre- 


dicted results of the current flight path. What was displayed 


to the MCO was the 


Current prediction of the impact point 


of the missile assuming that there would be no further control 


exercised over the 
Pitas ene center of 
mate of the target 
displayed relative 


is also rotated so 


missile. The display is generated such 
the screen corresponds to the last esti- 
position and the predicted impact point is 
to the center of the screen. The disvolay 


that motion of the impact point will 


correspond to control stick motion. Figure 7 illustrates 


the display screen 


system interface. 


and Appendix H further discusses the man/ 
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imi. LE SsiMULATION ORaRIE PROBEEM 


A proposed solution to a problem cannot be considered to 
be valid unless it is tested in some way in practice. One 
method of evaluation is to build an actual system and test 
it in a practical operating environment. This is generally 
too expensive to be feasible. Another method of testing a 
proposal is to simulate the system on an appropriate computer. 
The second method was chosen to evaluate the problem solution 
proposed in this paper. 

This simulation was done on two different computers. The 
first was an IBM System 360 Model 67, operating under the MVT 
Operating system. This system was used for initial program 
debugging since the Operating System has an excellent error- 
monitoring system. The second system was a hybrid computer 
consisting of an SDS 9300 digital computer and a COMCOR CI-= 


mO00 analog computer. 


A. SIMULATION OF TARGET LOCATION 
The simulation of the target location was accomplished 
by uSing the results of work in passive target location done 
by Demetry [l]. The results of the hybrid location method 
discussed in Chapter II were used by making use of the calcu- 
lated mean-square error values in target X and Y location 
as a function of time, from the hybrid location algorithm. 
Using these values involved a process of both normaliza- 


tion and interpolation. Since the simulation is designed for 
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a time increment of one second, the target location mean- 
square errors had to be interpolated since they were computed 
on a time increment of two seconds between .observations. 

Tne interpolation was made on the assumption that the 
mean-Square errors of location were piece-wise linear over 
the measurement interval involved. The normalization process 
was necesSary since when the values were generated, the air- 
craft was flying a set course which had a direct effect on 
tne value of the mean-square error. Thus, the values were 
normalized by multiplying them by a factor which was the 
closest point of approach (CPA) to the target for the prob- 
lem divided by the CPA to the target in the mean-Square-error 


CoMmpuracion program. 


Bee lL UGALLION OF TEE Jbl RE CONTROL SOLUTION 

A real-time simulation of the fire-control solution was 
eliminated after some research into the problem. A constant 
launch angle was used in its place. The solution was arrived 
at py assuming the missile to be a drag-free ballistic mass 
with the propulsion system of the true missile attached. 
Tnis approximation was used in an iterative method of comput- 
ing the angle which would give the maximum flight down range 
for a particular launch altitude. The iteration was performed 
on both the launch altitude and the launch angle. For this 
investigation the angular increment was taken as one degree 
and tne altitude increment was fifty meters. 

For each launch altitude, the angle which gave the maximum 


down-range flight of the assumed missile was plotted against 
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the launch altitude. The results were surprising in that it 
was found that for altitudes between fifty and five thousand 
meters, the angle varied only between 47 and 48 degrees. The 
launcn angle is defined as the initial angle between the hor- 
i1zontal and the velocity vector. This initial angle is re- 
duced by tne effect of gravity to an angle at burnout which 
will give the maximum down-range flight of the mass. The 
force of gravity is thus responsible for the initial angle 
being greater than 45 degrees. 

It must be noted that this angle was chosen to give the 
maximum down-range flight of the assumed missile. Thus, by 
shooting the true missile at this angle the maximum possible 
down-range flight would be possible for zero control. By 
the application of control to the missile in the pitch mode, 
this range can be either increased or decreased. In the 
desired situation this range should have to be decreased by 
using down pitch. This action would cause an increase in 
tne A/R ratio which has already been demonstrated to be 


desirable. 


©. SIMULATION OF @PHE MESSILE 

The missile model was the most important part of the simu- 
lation and thus was treated with the most rigor. The attempt 
was made to use a fairly accurate model in that dynamics 
approximating an actual missile were assumed. Appendix F 


contains a description of the missile model. 


a3 


Dee oUlDANGE “oaYNTHESIS 
The guidance synthesis for the simulation was done in two 
BeaGwes; one part for fhe simulation while i: “was execured@am 
tne System 360 and the second for the rt data-gathering 
Simulation on the SDS/COMCOR hybrid. 
1. System 360 Guidance Synthesis 
The guidance synthesis used during the preliminary 
runs on tne System 360 was a primitive, bang-bang control. 
Yaw was either plus or minus two degrees per second and pitch 
was a plus or minus 0.5 radian deflection of the control sur- 
face. The assumptions on which this control was synthesized 
were tne same as those discussed in the following sections. 
2. SDS/COMCOR Guidance Synthesis 
Before going into the discussion of the methods used 
to develop the guidance commands in the hybrid simulation, the 
assumptions on which the command synthesis was based will be 
discussed. 
as Rell-™raecpendent Comtrol 
The assumption that control is roll independent 
is valid since the majority Of missile control systems operate 
in this manner. Also the assumption allows a simplification 
of the command vector synthesis. The missile is thus consid- 
ered to either be roll stabilized or the control system in 
the missile is so designed to compensate for any roll of the 


missile. 
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bb. Two-Componéneseemerol Vector 

With roll independence assumed, the control vector 
consists of two components, a pitch and a yaw component. 
These components were limited in their magnitude by the syn- 
thesis program. 

The yaw was limited to a maximum rate of twenty 
degrees per second by the control mechanism itself. The 
pitch was limited to plus or minus forty-five degrees of 
control deflection. Additional limits on the control deflec- 
tions were built into the missile model. Yaw rate was lim- 
ited so that radial acceleration on the missile would not 
exceed ten G's. 

c. Control Vector Formation 

The control vector was formed by the human opera- 
tor directing a control stick. The maximum deflection of the 
stick corresponds to the maximum amount of control which 
could be generated. This control was then applied to the 
missile each second until the control stick was neutralized. 
The control affected the missile in that the control surfaces 
of the missile remained deflected only as long as the stick 


was held in a non-neutral position. 
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The simulation of this problem was performed only a few 
times with the completed program. As a result it was not 
possible to draw concrete conclusions about the validity of 
the total concept tested. However, several observations were 
made about the program and some areas where further investi- 


gation could be made were discovered. 


A. OBSERVATIONS 

One observation from the running of the program supported 
tne supposition that the use of impact vorediction in an air- 
to-surface environment was effective. The errors in impact 
prediction seemed to be related closely to errors in producing 
the measured missile state vector. Overall the impact pre- 
diction approach to missile control appeared very feasible. 

Another observation which was made concerned the initial 
launch angle assumed for the simulation. Apparently the use 
of a fixed angle of launch intended to give maximum down- 
range flight was not as good an.asSumption as first thought. 
For targets which were close to the launching aircraft, it 
was possible to overfly the target by shooting for maximum 
range. It would be much more desirable to have the launch 
angle be a function of the initial range to the target. 

The overfly or fly-by effect can be a serious problem. 
If the missile bypasses the target with a small angular offset 
to either side, it 1s possible that immediate application of 


yaw control in such circumstances will cause the missile to 
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circle the target at a constant range. This can occur since 
the missile is limited in the amount of radial acceleration 

it can withstand. It should be possible to overcome the prob- 
lem by the use of tactical doctrine if the missile does over- 
fly the target. One suggested doctrine would be to continue 
the fly-by course until the missile is far enough out to make 
the turn and hit the target area. This maneuver is illus- 


trated in Fig. 8. 


Pee PORTHER INVESTIGATION 

Three areas exist which the author feels merit further 
investigation. The areas are 1) prediction accuracy and 
methods; 2) radar filter accuracy; and 3) display methods. 

The accuracy of impact prediction should be investigated 
further. Although it was shown that the impact prediction 
method worked, no study of the accuracy of the method em- 
ployed was done. Also, it 1s possible that different meth- 
ods of impact prediction would prove to be more accurate 
than the velocity-proportional model used in the simulation. 

Prediction accuracy, from all indications, was very 
dependent on the accuracy with which the missile state was 
estimated. The Key percion Of tNe program In the missus. 
state estimation was the radar vector filter. The filter 
used in the problem was based on the assumption of a constant- 
velocity model of the missile. The use of a filter based on 
constant acceleration should be investigated as lags were 
noted in the measured missile state compared to the true mis- 
Sile vector when the missile was making hard or continued 


Maneuvers. 
SF 


The method of displaying the prediction information Ee 
the MCO also deserves further investigation. The display as 
1t was implemented did not provide any indication of the veloc-— 
ity of the missile. It is felt that if the MCO were given 
missile velocity information, he would be able to better con- 
trol the flight of the missile. This is especially true in 
the fly-by case. 

As the missile bypasses the target, the display gives a 
false indication as to the control necessary to hit the tar- 
get. This problem is generated as a result of the almost 
180-degree change in angle by which the display is rotated. 

Figure 9 illustrates an improved display which should 
solve the rotation problems caused by fly-by and also give 
velocity vector information to the MCO. The display is of 
the entire area of the problem. The axes of the new display 
are neither translated nor rotated as they were in the 
Simulation. 

Three points are included on the display. One is the 
current XY meaSured position of the missile. This is the 
left-most point in Fig. 9. The other two points are the 
current predicted impact point and the last estimate of the 
target location. A line from the missile position to the im- 
pact point will indicate the velocity information desired. 
Also the non-rotation of the display axes should alleviate 


the false indication of necessary command. 
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Missile Path 


L 


Figure 8. Course for Fly-by Error Correction 





Figure 9. Suggested Improved Display 


So 


om Roe ksy Guam ses scOlt 

The figures in this section show the results of only one 
run of the simulation program. The figures give two planar 
plots of@the fligmM:e path. One is a plot of the XY plane 
flight of the missile (Fig. 10) and the second is a vertical 
plane view of the flight path (Fig. 11). 

Other plots are included which show the control vector 
as a function of stime (Fig. 127 and the missile locatifon 
errors. as 4 LUnNCLiOn Ge timemtrig. isis 

It should be emphasized again that this is only one run 
and should not be considered to be in any way conclusive. To 
reach any definite conclusions a much more exhaustive investi- 


gation would have to be made. 
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APPENDIX A 


TARGET LOCATION BY SUCCESSIVE ANGULAR MEASUREMENTS 


FOR THE THREE-DIMENSIONAL CASE 


In chapter II there was given a brief description of 
locating a target in the two-dimensional plane by processing 
two sets of angular measurements and two associated aircraft 
positions. In this appendix the two-dimensional case will 
be expanded to the third dimension. 

Figure 15 shows the three-dimensional geometry involved 
in this problem. This figure is then shown in two plane views: 
the XY plane projection in Fig. 16 and the RZ plane projec- 
tion OF tR1g..17. 

From Fig. 16 the two locations of the launching aircraft 
are marked as xX, Yy and Xo, Yo with the associated azimuth 
angles 9, and 9.4. Using the geometric relationship for tan- 


a Z 


gent, the tangents of the angles shown are 


xX, 7X 

ae = tan Oy (A-L) 
eye 

X, 7X 

ee = tanh 04 e (A-2) 
ct 2 


Treating the equations as simultaneous equations and setting 
them up for olution, 


(A-3) 


i] 
~*~ 
! 
K 
ct 
Oo) 
2) 
D 


x, - Y, tan 4 l l 1 


Sees tan Ole xX (A-4) 


E t Z 2 2 2° 
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fetrang tan Oo) ee ea 04 = 8B, xy - Yy tan eT = C, and 


X, - Y, tan O5 = D and substituting these values into equa- 
tions A-3 and A-4 results in 


Ay = BY, =D. (A-6) 


Rewriting these equations in matrix form gives 





1 -aAl . t c | 
1 -B Y | . | — 
| 7 DI 
which has a solution of the form 
XI al ae Cj 
y “a 7 -B 'D (A-8) 
— | | : 


‘- = 
~ — 


The solution of this equation, arrived at by normal matrix 


algebra, 1s 


[xX i-B Al “| 
iE L 
= —— (A-9) 
Y 2 BI, 1 Di 
ic L 
Making substitutions back into the solution, the xy and Vy 
of the target are 
— tela 8, (X,-¥,tan6,) - tan 8. (X,-Y,tan6,) oe 
E tan 6. - tan @ 
a 2 
(X,-Y,tan6.) ~ (X,-¥,tan6,) 
fe tan 5 - tan 05 


The solution for locating the Z coordinate of the target 
will now be derived. 
Figure 17 will be used to illustrate the derivation of 


the target Z location problem. Again uSing the tangent 
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Figure 16. XY Plane View Location Geometry 
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Figure 17. RZ Plane View Location Geometry 
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relationships for the angles involved, it can be seen that 








ee b (A-12) 
D IL 
AL ; 
L~—-h 
2 ies igi : (A-13) 
D. Z 


From Fig. 16 the values for D, and D, can be seen to be 


1 Y 
WD Des 
Dy = ((X,-X, ) - (Y,-¥,) ) (A-14) 
a a Zn : ee _ 
ae ((X, X.) (Y, Y,) ) , (A-15) 


The fact 1s very evident from equations A-14 and A-15 
that once the X and Y coordinates of the target have been 
found, any single measurement of the depression angle may 
be used along with the Z location of the aircraft at the 
time of the depression angle measurement to determine the Z 


coordinate of the target. 
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APPENDIX B 


EFFECT OF A/R RATIO AND TARGET Z LOCATION ERROR 


ON RANGE MISS DISTANCE 


Several references have been made in this paper to the 
A/R ratio and its importance in determining missile miss dis- 
tances at the target. This appendix is devoted to a discus- 


sion of the importance of this effect. 


A. DEFINITION OF A/R RATIO AND ITS EFFECT ON MISS DISTANCE 

The A/R ratio, as defined for a parabola, is the ratio of 
the height of the parabola at its apex to the distance from 
a point under the apex to one of the terminal points of the 
Marabola.* This ratio is illustrated in Fig. 18. It Gan be 
seen from the figure that the larger the A/R ratio, the more 
nearly perpendicular the terminus of the parabola at its 
intersection with the axis. The more perpendicular the ter- 
minus of the parabola, the less the point of intersection 
with the target plane will be affected by an error in locating 
the target plane. Figures 19 and 20 illustrate the relation- 
ship just described. 

In Fig. LO ewriren.ilTUstraces™ a small a7 ReGak ©) mileee om 
me seen that for an error in starget 4 location 4, saeco. 
tively large error in range, AR, is generated. The alternate 
and desirable case of the large A/R ratio 1s shewm Inet ia. zoe 
For the same error in target Z location as in Fig. 19, the 
range error, AR, is much smaller than in the small A/R ratio 
example. 
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Figure 19. 


Figure 20. 
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B. CURVES RELATING MISS DISTANCE AND A/R RATIO AS A FUNCTION 
OF TARGET AZ. 


The computation of a family of curves relating miss dis- 
tances and A/R ratio as a function of target AZ 1s stralghiiee 
forward. Assuming a ballistic, drag-free mass, the time for 
that mass to fall the distance A under the influence of 


gravity was determined. This was done for the initial condi- 


tion of the mass having a zero Z velocity. In equation form 
t, = ees (B= 
b 
where Z = -9.8 Vee" and A is the initial Z of the mass with 


the target Z assumed to be zero. By using this time and the 
knowledge of the distance R from the apex of the parabola to 
the target, the initial R velocity was computed as 


Vo=— . (B-2) 


This velocity was then set as constant throughout the 
current run in which the Z location of the target was varied 


away from the perfect zero location for which the R velocity 


was computed. For each value of AZ, a time to fall, ter was 
computed and then substituted into the equation 
R=V te - (B-3) 


The difference between this distance thus computed and the R 
to the target was the miss distance. The series of computa- 
tions was repeated for a number of different A/R ratios and 
ave Si). 

The results were then plotted with A/R ratio on the X 


scale and miss distance on the Y scale. The family of curves 
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Ge 1G. 


21 was generated by connecting the miss distances at 


the different A/R ratios which were the result of the same 


nL. 


The values used in this problem were 


Rg = 


5000 meters (constant) 


A‘S such ‘that sA7R scab Smet 5 dig 2p pee ee 
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Figure 21. 


APPIN DT Xe 


DERIVATION OF THE KALMAN FILTER 


FOR THE MISSILE-TRACKING RADAR[2] 


In chapter II the Kalman filter for the system radar was 
discussed. This appendix gives the assumptions which were 
used in the derivation of the filter. 

fa eee was designed with the assumption that the 
missile-tracking radar provides measurements of range, range 
rate, azimuth angle, and depression angle to the missile 
relative to the aircraft. Another assumption was that the 
measurements were contaminated by additive noise whose sta- 
tistics were known. It was also assumed that the missile was 
not subjected to any random perturbations. Based on these 


assumptions, some of the matrices of equations C-l, C-3, and 


C-5 may be defined. 


ey, = Py jkngtt Mey ip (C-1) 
Eines (T-G)H) Py Ray (C22) 
Pelt k = OPE Ke + Q (eae) 
Xe | = eas Gye (Zy HX. yy) C4 
Xelker = Ser feed (C-5)* 


The Q matrices are identically zero since it was assumed 


that the missile was subject to no random perturbations. The 


the symbols used in the equations C-l through C-5 and 
elsewhere in this appendix differ somewhat from the svmbol- 
ogy used in R. C. K. Lee. [Ref. 2]. 


Suh 


R matrices are not zero since they convey to the filter infor- 
mation regarding the variance of the additive noise in the 
measured values. 

The ¢ matrices are the same for all the filters since the 
same dynamics (constant velocity) were assumed for all the 
filters. Equation C-6 gives these ¢ matrices as 


ioe 8 
® = = == : (C-6) 
O dL; 


The H matrices in equations C-l, C-2, and C-4 are measure- 
ment matrices which determine which variables are measurable 
by the filters. In the range filter HG is a2xX 2 identity 
matrix. The Hs and ae are 1 X 2 matrices with the (1,1) 
elements being unity. 

It should be noted that equations C-l, C-2, and C-3 are 
data-independent as long as the matrices involved in computing 
them are data-independent and if the H, R, 9%, and Q matrices 
are constants. If this is true, then the gain schedule can 
be computed offline from the actual simulation of the problem 
Or operational environment and simply inserted into the com- 
puter memory for use at the proper time. With the precomputed 
gain schedule in use, the implementation of the filter to 
process data in real-time is a relatively easy matter and 
does not contribute greatly to the computational load on the 
GOmou te 29 

An important consideration in designing any filter is the 
expected value of the measured parameters which will be used 


to initialize the filter. This is especially important when 


choosing values for the initial covariance matrices, ac 
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As these filters were to be employed in conjunction with a 
missile-tracking radar, the assumptions were made that the 
initial range to the missile and range rate of the missile 
would be unknown and that the initial angles to the missile 
would also be unknown. These assumptions were reflected in 
the covariance matrix initializations when the gain schedule 
was computed. The values used are given for the range and 
range rate filter in equation C-7 and the angular filters in 


equation C-8. 


10000. m*@ 0 | 

P = (c=?) 
L}0 0 1000. m“/sec* 
10 rad- 0 

P = | (C-8) 
1}0 0 iOMEeNCIg eee | 


The values used for the R matrices in the gain schedule 


computation were 


100 ne 0 
= (ooo) 
rg 0 ino m?/sect 
R= R, = .0004 rad? (C-10) 
az dp 


By using very large initial covariance matrices and seu 
ting the initial estimates of the state vectors to zero, the 
gain element (1,1) for the first time increment will be very 
close to unity, thus causing the first measurement of the 
particular state variable to become the first estimate of 
that variable. As time passes, however, the incoming meas- 
urement is weighted less and less as the gain goes down, 1.e., 
the filter becomes more confident in the accuracy of its 


estimates. 
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However, this reduction of gain as a function of time can 
ities be self-defeating. In the case in which this filter 
was employed, the reduction of the gain with time could be 
disastrous. Some corrective action had to be taken to prevent 
the gains from approaching zero and, as a result, having the 
filter ignore new data. Otherwise, as the missile maneuvered, 
the filters would have very quickly begun to produce very bad 
estimates of the radar state vector. 

Thus, the effect which is deSired is one which will not 
allow the gains to approach zero. There are two actions which 
will prevent this from occurring: letting the Q matrices be 
non-zero, Or truncating the gain schedule after some length 
of time kT. Since the Q matrices have already been considered 
to be zero, there remains gain-schedule truncation. By trun- 
cating the gain schedule the filters continue to accept in- 
coming measurements as containing significant information, 
and the amount of core storage in the computer is reduced 
greatly. 

The gain schedules in this particular simulation were 
truncated in such a manner that the previous estimate produced 
by the filters were weighted by a factor of approximately 0.8 
and the incoming measurements were weighted by a factor of 


approximately 0.2. 
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APPENDIX D 


DERIVATION OF THE SOLUTION Leer 


IMPACT PREDICTION EQUATIONS 


This appendix discusses the derivation of the impact 
prediction equations used for the missile impact point pre- 
diction. Two methods are discussed: one in which the drag 
force 1S proportional to the velocity and a second in which 
the drag force is proportional to the velocity squared. The 
first was used in the simulation program. 

In order to derive closed-form solutions to the equa- 
tions of motion for the missile, it was assumed that the drag 
coefficient was constant over the prediction interval. 

This assumption will introduce errors into the predictions. 
However, this error should be reduced the closer the missile 
approaches to its terminal point. Near the terminal point of 
the flight the greatest error will be in the predicted time 
until impact since the Z acceleration will be the greatest 
at that time. This should not greatly affect the X and Y 
prediction since the X and Y accelerations will tend to be 
very small at that point. 

The method used to solve the equations is to first solve 
the Z motion equation for the time at which the missile Z is 
equal to the target Z. The value of time thus found, which 
is the predicted remaining flight time of the missile, is 
then used with the current X and Y positions and velocities 


to compute the predicted impact point. 


OnE 


hae METHOD IT — DRAG PROPOREEGNAL TO VEE@GrETy 

In method I the drag force on the missile is considered 
to be proportional to the velocity of the missile. The equa- 
tions for this case are linear differential equations. 


The equations of motion for method I are 


mX = -K,X (Do 
mY = -K5Y )) 
mz = -K,2 - mG (D-3) 
where 
K, = o[ XA (D-4) 
K, = o|¥| A, (D-5) 
K, = 0/2) A (D-6) 


in which the x, a and Z are the current values of the 
velocities of the missile, p is the atmospheric density, and 
A. is a reference area of .114 meter’. 

The solution for the Z motion equation D-3 has the form 


Sie 


0 = -Ae -Bt+ec (D~7) 
where 
Tig t 0) 2 
Az Mm 4 ms (D-8) 
Ne K 
3 
B= ue (D-9) 
3 
C= A+ Z (0) - Z (D-10) 
K 
a . (ese 
m 


The problem now becomes the solution of the time-domain 
equation of the Z motion for the time at which the missile 2 
is equal to the target Z. The C coefficient in equation D-10 
is set to reflect this desired equality already. 

The method chosen to solve this equation is the Newton- 
Raphson method. Upon commencement of the iterations all of 
the required coefficients were computed using the current 
values of the positions and velocities. The termination test 
constant used to terminate the iterative process was 0.1 sec- 
ond. Whenever one value of time had changed from the previous 
value of time by an amount less than or equal to the test 
increment, the iterative process was terminated. The value 
of 0.1 seconds was chosen for the test value since regardless 
of the sample interval given to the simulation program the 
missile model was advanced in increments of 0.1 seconds. 
Thus, whenever the time was less than or equal to the minimum 
increment that the missile could be advanced, the iteration 
was stopped. 

The value of time which was returned from the Newton- 


Raphson iteration was then substituted into the equations for 


X and Y impact point prediction. These equations are 
mx (0) i 
coy. . Ce ge oa) 
p m Ry 
mY (0) -vt 
es = Y 69) a K, (1.0-e ) (D-13) 
where 
K 
m 


od 


al 
NO 


(D=-15) 


These are the equations for method I prediction. The 
method II equations are more complicated to solve since they 


are non-linear. 


B. METHOD IIT - DRAG PROPORTIONAL TO VELOCITY SQUARED 

In this derivation of impact prediction equations the 
drag force on the missile is taken to be proportional to the 
Square of the missile velocity. Again, the drag coefficients 
are computed at the beginning of the solution and considered 
constant throughout. 


The equations of motion are 


mX = -K,X (D-16) 
oo e p 
mY = “K5Y (D-17) 
ee @ 2 
mz = “K 4 - mG : (D-18) 


The steps of this solution are the same as used in the 
previous method: 1) solve the Z motion equation for Z asa 
function of time, 2) solve for the time when missile and 
target Z are equal, and 3) use the time thus found to compute 
— and ue from the prediction equations. 


First solve for Z as a function of time, writing equation 


D-18 as 
avout 2 = 
where 
K 
C= _3 (D-20) 
m 
and 
y=2 : (D-21) 
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Separating variables and integrating, 


a Sane ee) te Dy =t. (D-22) 
¥CG 
At t = 0, y = ZO). There fore 
iE -l,, 
OT oe col (Z_ (0) ¥C/G) 3 (D-23) 
VCG 
Eemving for y and substi tuteng for Dy, 
y - YG/C tan(D,vCG) = -/G/C tan(tvCG) (D-24) 
Oi 
y = =) 6/7 CG, tamger CG) at Z (0) : (B= 25) 
Substituting y = ae and letting Z_ (0) be constant, 
dz = [-VG/C tan(tVCG) + Z_ (0) Jdt (D-26) 
and 
Z = -VG/C —— log, cos(tv6G) + 2 (0)t + Dz (D-27) 
/CG 
wee = 0, 2= Z (0). Therefore D. = 4)(0) and 
m 2 m 
_ A 2 
Z(t) = a log, efor (ere) sr Zi (O)t + Z (0) (D-28) 


where zZ_(0) and Z_60) refer to the Z velocity and position 
of the missile at the time the prediction is made. 

Solving for the time when Z(t) is equal to the target 
elevation is done by substituting Z(t) = Zo and using the 
Newton-Raphson method. 

The solution of the X and Y motion equations for impact 
point will be shown for the X equation only and the Y equa- 
tion will be written by inspection. 

ae 


mX = -K,Xx (Deo) 


ae (D-30) 
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. Separating variables and integrating 


Cy + Dy =t . (D=31) 


hea = 0, y= xX (0). Therefore 
it 








Ds = oe ee ° (D-32) 
CrG(.0) 
m 
ax ; 
SUubSsecrruring, y = as and letting X_ (0) be constant 
Cine _ : 
ia oe CD, (D— 35 
Separating variables and integrating 
xX = a vege cl = Siam + D : (D-34) 
C e 2 
X_ (0) 
m 
St j=) 0 ee ae xX (0). Therefore 
Deen (ON = aliens (= =) (D-35) 
2 m e Je 
X_ (0) 
m 
and 
X(t) = * log_(Ct - ) - = log_(- -< z ) + X_(0) (D-36) 
e c e m 
X_ (0) ane!) 
m m 
i 1 ; 
X(t) = & log (Ct - a——) (-X_(0)) + X_ (0) (D-37) 
& m m 
X_ (0) 
Or @ 
ze K, tx (0) 
hee) = zo log, [1 ~ ]} + X_ (0) : (D~ 38) 
For this method the prediction equations are 
8 K,tX_ (0) 
AS = K, log, [1 ~ ) + X_(0) (D-39) 
and ; 
a K,tY_ (0) 
* = ey log {1 = ] - 69) . (D-40) 
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APPENDIX E 
REFERENCE SYSTEMS FOR THE PROBLEM 


The reference systems for the problem are based on the 
position of the launch point of the aircraft involved. The 
inertial navigator of the aircraft is initialized when the 
aircraft is launched. All other points in the problem are 
also referred to this initialization point of the aircraft 
even though they may not be measured with respect to that 
point at first. The measured missile state vector is an 
example of this type in that it is first referred to the 
mother aircraft and then converted to a vector which is 
referred to the initialization point. 

The cartesian reference is a standard three-dimensional, 
right-hand reference system. Whenever a reference is shifted, 
the new frame is also a three-dimensional, right-hand system. 
Figure 22 illustrates the various reference points and systems. 

Azimuth angle is defined as the angle between true north, 
positive clockwise to the line-of-sight from the aircraft to 
the item of interest. In order to define the depression 
angle, it is first necessary to define the normal plane. The 
normal plane is defined as being a plane normal to a line pas- 
Sing from the aircraft to the center of the earth. Depression 
angle is defined as the angle between the normal plane and the 
line-of-sight. This LOS must be contained in a plane which is 
perpendicular to the normal plane. These angles are shown in 


Fig. 23 and 24 respectively. 
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Figure 22. Problem Cartesian Reference Frame 
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Figure 23. Definition of Figure 24. Definition of 
Azimuth Angle Depression Angle 
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APPENDIX F 
PHS MISsticeh MODEL 


The missile model was relatively honest in its aerodynamic 
performance although it was not a true missile. The maximum 
yaw rate of the missile was limited to allow radial accelera- 
tion loads not to exceed ten G's. This was thought to be 
representative of current capabilities. 

The missile was powered by a rocket motor of 3000-lbs. 
thrust with a burn time of 30 seconds. After the burnout the 
missile glided for the remainder of its flight. 

The model was so designed that the mass and moment of 
inertia changed as the rocket motor burned. The total pro- 
pellant weight was approximately 163.3 kg. The reduced mass 
and moment of inertia were then used for the remainder of the 
flight. A free body diagram of the missile appears in Fig. 


25. The equations of motion are 


mR = -D cosy - L siny + T cosé@ (F-1) 
mz = -D siny + L cosy + T sin8 - mG (F-2) 
I 9 = -jL - 3D sino + 2CF - K 6 (F-3) 

VY J J 42 

where 

Onn ov" (F-4) 
D = Kj (F-5) 
L = K,Qa ee 
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and 


ic 
ey 


(thrust) 

flee, 

(drag) 

(veloe? ty) 

(Cont no erOrce) 
(velocity vector angle) 
(angle of attack) 
(missile axis angle) 
(axis rate) 

(axis acceleration) 


(mass) 


(gravity) 
(control moment arm) 


(center of pressure 
moment arm) 


(moment of inertia) 


(air density) 


(drag constant) 


(control force 
COnSstane) 


(damping constant) 


(lv it cometary 


13344. N 


m/sec 


rad 

rad 

rad 
rad/sec 
2 


rad/sec 


453.6 kg 
(JIU. SGmetcer Olt nem) 


O26 m/sec? 


Zi 


Loom 


Blo youre ia 
(ON 7 eat ae Ceactaatenteneey 


1.284 kg/m? 


00556 nt/(kg/m>) (m/sec) * 


C0 nt/ (kg/m) (m/sec) “rad 


,004 nt/(kg/m?) (m/sec) * (rad/sec) 


.008 nt/(kg/m>) (m/sec) “rad. 
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APPENDIX G 


THE LAUNCHING AIRCRAFT MODEL 


The launching aircraft model was a very simple model for 
the simulation. The aircraft was given an initial heading 
and velocity. It was then assumed to fly a straight line 
making no turns. The location of the aircraft was updated 
each sample period. 

Although preliminary investigation had shown that the 
most critical component of noise in the entire system was 
that in the inertial navigation package, the inertial navi- 
gator was considered to be noiseless. The reasons for 
neglecting the noise were that 1) the primary purpose of this 
investigation was in the field of application of impact pre- 
diction to an operator-controlled, air-to-surface weapon 
system, and 2) no figures were available concerning sample- 


to-sample noise in inertial navigators. 
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APPENDIX H 


MAN/SYSTEM INTERFACE 


In any system with a human operator, it 1s important that 
information be presented concisely and in such a manner that 
it can best be used by the operator. The communication link 
in any such application is a two-way channel. The system 
sends information to the operator in one direction and the 
Becrator-generated control commands, trave Min) Enewe ego: 

the link to tMe operator Carried two items of informmavies: 
The first was the most recent prediction of impact point and 
the second was the predicted time until impact. 

The first item, impact point prediction, was displayed 
to the operator on a display scope with the target assumed to 
be in the center. This presentation was limited in such a 
manner that the predicted impact point could never appear 
off-scale. Figures 26, 27, and 28 give examples of the 
limiting action which was employed. The display, when in a 
Limited state, gave no information as to the error of the 
predicted impact point except the direction in which it was 
off-scale. 

The display was generated such that the missile always 
approached the target from the bottom of the scope. MThis 
required a transformation of the axes of the display. This 
was necessary so that the motion of the control stick always 


drove the missile toward the target. 


12 





Figure 26, Impact Point Within Scale 
on Display 





Figure 27. Impact Point Off-Scale 
Left on Display 





Figure 28. Impact Point Off-Scale 
Left, Off-Scale Vertically 
on Display 
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The transformation of axes consisted of both a translation 
and rotation. The translation was accomplished by generating 
a reference frame which had its origin located on the target. 
The predicted impact point was then referred to this origin 
and the axes rotated in such a manner as to make the Y axis 
of the rotated system coincide with the LOS from the current 
measured missile position to the target. This process is 
meliustrated imehig. 29 and 20n 

It can be shown that the location of the predicted im- 
pact point in the translated and rotated reference frame 


(display reference system) is 


= ne (H-1) 
where 
x 
a a 
Xa = ; (H-2) 
d 
cos 8 -sin 8 
Sin 6 cos 98 
X - X 
= P t = 
*d¢f (H-4) 
Y - Y 
p iE 
and 
= cae J . = 
6 = Seam C(x, X SY, Yi)? : (H-5) 


The second item of information presented to the operator 
was the predicted time remaining until impact. This value 
was computed once each sample period throughout the missile 
flight. It was generated on the digital computer, converted 
to an analog voltage, and displayed on a digital voltmeter. 


The display was scaled such that times were displayed only if 


J 
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Figure 29. Translation of Axes for Display 





|X 


Figure 30. Rotation of Axes for Display 
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they were less than or equal to 100 seconds. Any value of 
time which was greater than 100 seconds was displayed as 100 
seconds. 

Commands from the operator to the system formed the other 
direction of the two-way communication link. These commands 
were generated by a control stick which the operator used to 
guide the missile. Full deflections on the stick corresponded 
to full deflections on the missile controls. The pitch con- 
trol was scaled to provide 45 degrees of control movement on 
the missile when the stick was fully deflected. The yaw con- 
trol provided a yaw rate of 20 degrees/second when fully de- 
flected. This value was also subject to limiting by the 


missile dynamics. 


ca 


APPENDIX I 
AEDES CUDOSZON OF IHL SIMUBATION PReerRaM 


Tnis appendix describing the program used in the simula- 
tion of this problem is included in the event that some other 
investigator would desire to continue research in this par- 
ticular field. In it the author has tried to include all 
peiGtinent inmormablem concerning the routines and their )fiies 
tions and especially the programming "tricks" which were used 
in some instances to provide a particular function in the 
overall program. 

The appendix is divided into four sections. The first 
section 1S a summary of the routines and the function which 
each performs. The second section lists the routines and 
the variable names within each routine. The method of inter- 
routine communication is the subject of the third section 
and the programming "tricks" used are discussed in a fourth 
section. A complete listing of the program and an analog 


patching diagram follow the appendix. 


A. SUMMARY OF THE ROUTINES AND THEIR FUNCTIONS 
1. Main Line Program 
The main line program serves to perform two functions. 
The first is as a vehicle for input to the program and output 
From the program, and the second is as a calling and control- 
ling sequence for the various routines which do the actual 


computation for the simulation. If operated in a segmented 


v3 


mode of core allocation, the main line program is always core 
resident. 

several test functions are performed by the main line 
program also. It testswfor the condititon of preblen@idemi— 
nation. Problem termination occurs when the measured Z of 
tne missile is equal to or less than the target Z. It also 
controls switching from the passive target location mode of 
the problem to the active mode in which the missile is fired 
and guided to the target. 

2. Initialization Routines 
a. Inertial Navigator Initialization 
The subroutine INIT initializes the inertial 

navigation routine by setting the initial position of the 
aircraft including the altitude. It also initializes the 
heading and computes component velocities. The final items 
initialized by the routine are the elapsed time counter, 
missile flight time counter, and the problem time increment. 

b. Missile Initialization 

The subroutine MINIT initializes the missile model 

at the time of launch. All parameters of the missile model 
are set equal to the aircraft parameters at the time of launch. 
Also the routine gives initial values to the measured missile 
vector to enable other routines to function properly during 
the. frrst second@et the sami atiowsiirri1 Mall roWeimes Fane 


pene rating full outpurs. 


wo 


Pe ALrCraGevot mula ronmpRomueane 
The subroutine NAVPCK is the routine which provides 
aircraft state vector information for the problem. The 
routine provides a straight-line flight path for the aircraft 
at any initial heading from 0 to 359 degrees at any initial 
velocity. The routine receives flight path information in 
R/8 form and provides flight-path information in the three 
cartesian coordinates. 
4, Passive Target Location 
The subroutine PFIL simulates the action of the hybrid 
filter presented in section II-A-3. This location action is 
simulated by interpolating and normalizing values of variance 
provided for the hybrid method and generating noise which is 
added to the true target location. The final noisy location 
of the target is then used throughout the remainder of the 
prop lems 
5. Missile Model 
The missile model routine is of sufficient importance 
that an entire previous appendix was devoted to it. The reader 
is referred to Appendix F for a discussion of the missile model. 
6. Routines Involved in Missile Tracking 
a. Radar Vector Generation 
The subroutine REVCON is the routine which gener- 
ates the true radar vector from the geometrical relationships 
between the true missile and the aircraft. This vector rev- 
resents the true radar state before any noise is added and 


before quantization takes place. 


80 


b. Noise Addition and Quantization of Radar Vector 
The subroutine RADAR2 takes the true radar vector 
produced by REVCON and performs two operations on it. First 
the routine adds the measurement noise associated with the 
particular measured state variable. Second, the routine 
quantizes the measurements into the discrete form which 
would be provided by a digitized radar output unit. 
c. Radar Vector Filtering 
The subroutine AFIL provides the filtering of 
tne noisy, quantized radar vector to produce a radar vector 
whose variance is less than that of the noisy vector. The 
subroutine uses a schedule of precomputed filter gains and 
a switch which determinés the number of samples after which 
the gain schedule is truftated. This value can be from one 
to thirty seconds assuming gains computed for a sample time 
Strmeme Second. 
d. Spherical t6 Cartesian Coordinate Conversion 
The subroutine CORDCN processes the filtered 
radar vector spherical values to produce a measured missile 
vector composed of missile position and velocities relative 
to the true problem ref@rence. The subroutine first pro- 
duces a cartesian vector relative to the aircraft. The air- 
craft component positions and velocities are then added to 
the relative components to produce the final measured missile 


vector. 
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7; piimtial Missale GContwel d«oukine 
The subroutine SHOOT controls the flight of the mis- 
sile during the boosted phase (first thirty seconds) of its 
flaght. Theseroutine attempts to drave the missile velocity 
vector to assume an up angle of 47.5 degrees and to turn the 
missile velocity vector projection on the XY plane to coin- 
cide with the LOS from the missile to the target. 
o Sinpact Predvetion andebesp lay, 
a. Impact Point Prediction 
The subroutine PREDCT develops the X and Y im- 
pact point prediction values. This routine employs the sub- 
routine NR to assist in this task. 
b. Time-Until-Impact Solution 
The subroutine NR solves for the time to go 
until impact for the missile by using the Newton-Raphson 
iterative method of solution for differential equations. 
c. Command Generation 
The subroutine COMGEN converts the outputs of 
the control stick to values which can drive the missile mod- 
el. This routine uses the SDS 9300 routines for analog-to- 
digital conversion. Full-scale values of pitch correspond 
to plus or manus 45) degrees of control deflection and full- 
scale values of yaw correspond to turn rates of plus or 
minus 20 degrees/second. 
pee” KeGac tl Onmpisp lay manduStabrstuesmeROMEInes 
a. Display Generation 
The subroutine DISGEN converts the predicted im- 


pact points into values which can be used to drive the 


Bie 


Situation display. 


and translation necessary to generate a useful display. 


This routine performs the axis rotation 


The 


limiting of the display range is also performed. 


Jone 


Minimal Statistics Routine 


The subroutine MISS provides minimal statistics 


about the success of the run at the time of flight termina- 


1 ON « 


impact time of the missile, 


and total miss distance. 


This routine provides the X and Y miss distances, 


More 


complete statistics are generated by a separate program which 


analyzes data which is written on a storage medium during the 


execution of the program. 


The separate analysis routine pro- 


vides such results as missile flight profile, control applied 


vs. time, missile location errors vs. 


time, and others. 


B. ROUTINE SYMBOLS AND THEIR MEANINGS 


i. 


Main Line Program 
SyyAmlorsil 


GBQ/KBO 

RUNS /RNS 

NP 

ADP 


SMU, ARNG, ARNGD, AAZ, 


VL, THI, Z, TME 


ETIME, TIME, FTIME 


ETA 
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Meaning 


Time for active filter 
Galina euncation. 


Number of runs for the 
program. 


SGwltch uSeC@GmbenmentelLallze 
AFIL. 


Mean of noise and noise 
variances. 


Initial velocity, heading 
(in degrees), aircraft 
altitude, and time incre- 
mete 


Elapsed time, time incre- 
ment, flight time. 


Initial launch angle. 


Subroutine INIT 
Symbol 


XA,YA,ZA,XDA,YDA,ZDA 


Calne ed tN 


Subroutine NAVPCK 


Symbol 


ok Gab 


XA,YA,ZA,XDA,YDA,ZDA,TH,V 


Subroutine MINIT 
Symbol 
YAW , PITCH 


XMM ,YMM,ZMM,XDMM,YDMM, 
Z,DMM 


XM,YM,ZM,XDM,YDM,ZDM, 
XDDM,YDDM,ZDDM, 
HEAD ,GAMMA,THET, 
THETD,THETDD 

Subroutine PFIL 

Symbol 

Day. 


OnenOM., KM wd 


XS,YS,A 


MT VT CA 
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Meaning 

AlrGémaftt. location wang 
velocity in cartesian 
coordinate frame. 
Conversion factor (degrees 
to radians), converted 
heading angle, heading in 


degrees, velocity in 
meters/second. 


Meaning 

Conversion factor (degrees 
to radians), converted 
heading. 


Aircraft state vector. 


Meaning 
Missile feommandvectom 


Measured missile vector 
in position and velocity. 


True missile vector. 


Meaning 
True target location. 


Constants used for 
interpolation. 


X variance, Y variance, 
noise value returned by 
GAUSS. 


Noisy target location, 
CPA in problem. 


Subroutine MISAD 
Symbol 


RHO,T,G 


REFAR 


K1,K3,K4,K5 


K,YINC,ANG, VEL, VELSQ 


DELTA, GAMMA, ALPHA 


D,QL,QM 

QMASS ,QLYY ,QMULT 
Subroutine REVCON 
Symbol 


ae 


RPRI ,RDPRI ,T1 ,T2,VLXY 


R,RD,T,P 
Subroutine RADAR2 
Symbol 


el 


ITEMP , AN 


R,RD,T,P 
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Meaning 


Air density, thrust, 
epetshi Lie? 


Aerodynamic reference 
area. 


Drag, lift and moment 
constants. 


Advance count, yaw incre- 
ment, current missile 
heading, current velocity, 
current velocity squared. 


Control angle, missile 
axis angle, angle of 
attack. 


Drag, 11ft, and moment 
forces. 


Missile mass, moment of 


inertia, constant relating 
mass and moment. 


Meaning 
Work variables. 


Work variables for develop- 
ing the radar vector. 


True radar vector. 


Meaning 


Bin size for 13-bit shaft 
encoder. 


Temporary variable used 
in the conversions, and 
noise value returned by 
GAUSS. 


Noisy, quantized radar 
VeCEOr. 


ie 


dels 


Subroutine AFIL 
Symbol 


NE 


ZR1,ZR2,ZAl,ZD1 


XRI -XR2 pee D 1 , XD2 


ERG CAZAGDE 


Subroutine CORDCN 
Symbol 


CSP,SNP,CST,SNT 

CSPSNT,CSPCST,SNPSNT, 
SNPCST 

R,RD,T,TD,P,PD 

XMM ,YMM, ZMM 


XDMM,YDMM, ZDMM 


Subroutine SHOOT 
Symbol 


GAMMA, HEAD 


ETA,AL 


Mar, PL PCH 
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Meaning 


Switch to comtsol fileeem 
mya plea WicZ a EeeOln 


Measurements of range, 
range rate, azimuth, and 
depression angles. 


Estimated values of range, 
range rate, azimuth, 
azimuth rate, depression, 
and depression rate. 


Range,azimuth, and de- 
pression filter gains. 


Meaning 


Sine and cosine of 9 
ac, Gi 


Cross products of sine 
and cosine values. 


Range, range rate, theta, 
theta rate, phi, and ohi 
mouee 


Measured missile cartesian 
position. 


Measured missile cartesian 
velocity. 


Meaning 


Current missile velocity 
vector elevation and cur- 
rent missile heading. 


Desired missile velocity 
vector elevation and cur- 
rent LOS angle from the 
missile to the target. 


Missile command vector. 


ili, 


ies 


eds 


ss 


Subroutine PREDCT 
Symoou 
QMASS ,RHO, REFAR 


era 


XP,YP,ZP 


Subroutine NR 


Symbol 


QMASS,RHO,REFAR,G,V 


Kl ,ALPHA 


A,B,C,TOLD,TNEW 


lee ee 


Subroutine COMGEN 
Symbol 

YAW ,PITCH 

Numeric arguments 
Subroutine DISGEN 
Symbol 


XDF ,YDF 


XTA,YTA 


ALPH 
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Meaning 

Missile mass, air density, 
missile aerodynamic refer- 
ence area. 


Drag coefficients and 
missile velocity. 


X and Y impact position 


predictions, and time 
Until anbacu: 


Meaning 

Missile mass, air density, 
misSile aerodynamic refer- 
ence area, gravity, and 
missile velocity. 


Drag seoec f firemen basa 
exponential argument. 


Coefficients of prediction 
equations, last time esti- 
mate, new time estimate. 


Function value, and 
derivative. 


Meaning 
Missile Command vector. 


G@enVers tone facto cc 


Meaning 


Translated predicted 
IMpaGre eo limer 


Used to determine 
rotation angle. 


Rotation angle. 


MD yy D2.) Display points and time 
Unt eimoacc.. 


16. Subroutine MISS 

Symbol Meaning 

ZUM bly eT 2 Last preimpact missile 2 
position, times which 
bracket impact time. 

XIM1,YIM1,XI,YI Last preimpact X andi 
positions, and missile 
impact positions in X 
and Y. 


Ap, Een Impact time as elavsed 
time and flight time. 


Kl ,K2,QK2 Variables used to edit 
inaccuracies out of 
Piteie cime Clocw« 

MSM ST Wade X and Y miss distances, 


and last estimated target 
pes teElon: 


See PE RROUTING COMMUNICATION 

Inteérroutine communication is carried out throughout (tae 
program by means of labeled FORTRAN COMMON blocks. These 
blocks are designed in such a manner that variables relatina 
closely to each other have been grouped into sections named 
to reflect their functions. The block names are given below 
along with the variables included and a description of the 
mIMcGLon OL the bileek, 

1. Blank COMMON 

The variables included in blank COMMON are IX, ETIME, 

Pei ola | TMi was Originally s@tne Seed yalue for the 


GAUSS routine which produces random numbers. It was not re- 


tained when a change in random number generators was made. 
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ETIME, TIME, and FTIME are the problem elapsed time, time 
increment, and missile flight time respectively. 
2. Block TGPRD 
The variables included in the block TGPRD are the 
cartesian location of the target and the X and Y predicted 
impact points of the missile. 2ZP is the time until impact. 
3. Block SIGMAS 
This block contains the values of the sSigmas to be 
used in the radar to introduce noise in the measurements. 
These values are set in the main line program. 
4, Block GAINS 
The computed filter gains used by the active radar 
filter are stored in this dimensioned block. The block 
provides storage for 120 range filter gains and 60 azimuth 
and depression filter gains each. 
See Oo LOek bo LM 
This block contains the values of measured missile 
position and velocity developed by the missile tracking 
routine. 
Gr LOC mma c RET 
The values of aircraft position and velocity at each 
time k are held in this block. Also this block contains the 
aircraft velocity in meters per second and aircraft heading 
in degrees. 
(gee LOCK Missin 
This vector contains all values which pertain to the 
true missile. The vector contains position, velocity, and 
acceleration in cartesian coordinates; missile axis angle, 
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angular velocity, and angular acceleration; and velocity vec- 
tor heading and elevation angle. This is updated each sample 
period to reflect the true state of the missile. 
Bice Db EOC OMNI 
The commands to the missile in yaw and pitch comprise 
the two elements of this block. The yaw values are up to 
twenty degrees per second in either direction and pitch val- 
ues range from plus to minus forty-five degrees in control 
deflection. 
o eB Lock SrGs 
The values of variance which are used in the passive 
location routine PFIL are transferred to the routine by this 
block. Each of the values is dimensioned to accept fifty- 
nine values of variance. 
10. Block RADAR 
This block is used to transfer the values of the radar 
information through all of the routines of the missile track- 
ing portion of the program. The six elements of the vector 
are range, range rate, azimuth, azimuth rate, depression, and 


depression rate. 


D. PROGRAMMING "TRICKS" 

Only one programming irregularity exists in the simula- 
tion program. In the subroutine PFIL the variable ZP in the 
COMMON block TGPRD is used as a switch to allow the routine 
to hold the true target location after the first call. The 
variable ZP was chosen since it was already in the COMMON 
block and since it is not used by any other routine in the 


passive phase of the program. 
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